Sensitivity of Hawking radiation to superluminal dispersion relations 
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We analyze the Hawking radiation process due to collapsing configurations in the presence of 
superluminal modifications of the dispersion relation. With such superluminal dispersion relations, 
the horizon effectively becomes a frequency-dependent concept. In particular, at every moment of 
the collapse, there is a critical frequency above which no horizon is experienced. We show that, as 
a consequence, the late-time radiation suffers strong modifications, both quantitative and qualita- 
tive, compared to the standard Hawking picture. Concretely, we show that the radiation spectrum 
becomes dependent on the measuring time, on the surface gravities associated with different fre- 
quencies, and on the critical frequency. Even if the critical frequency is well above the Planck scale, 
important modifications still show up. 
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I. INTRODUCTION 

Phenomenological approaches to quantum gravity have 
recently started to develop in parallel with the more tra- 
ditional attempts to construct such a theory from first 
principles. In particular, increasing attention has focused 
on the consideration that maybe Lorentz invariance is not 
a fundamental law, but an effective low-energy symme- 
try which is broken at high energies (see, e.g., [1] for a 
general introduction). Conceptually speaking, in quan- 
tum gravity theories from first principles, it is not re- 
ally clear whether Lorentz invariance is fundamental or 
effective, and in the latter case, how its breaking scale 
is related to the Planck scale. For example, while many 
string theory scenarios axiomatically incorporate Lorentz 
invariance, it has been argued that in certain situations, 
violations of Lorentz invariance may occur in a way con- 
sistent with world-sheet conformal invariance 0], thus 
leading to acceptable string theory backgrounds. In the 
context of loop quantum gravity, in Q it has been ar- 
gued that quantum effects should modify the relativistic 
dispersion relations, although the issue seems far from 
settled (see for some general remarks). In scenarios 
of emergent gravity based on condensed matter analo- 
gies Q, the situation is clearer: Lorentz invariance is 
a low-energy effective symmetry, and so it is expected 
to break at some scale, although not necessarily related 
to (and therefore possibly much higher than) the Planck 
scale 0- 

A simple way of modelling a wide range of Lorentz 
violating effects (and a quite natural one in the case of 
condensed matter analogies) consists in modifying the 
dispersion relations at high energy This modifica- 
tion can be subluminal or superluminal, depending on 
whether high-frequency modes move slower or faster than 
low-energy ones. We should mention that there also ex- 
ist ways of modifying the dispersion relations, for ex- 
ample by introducing a minimum length or maximum 
energy, without violating Lorentz invariance Q. Here, 



however, we will be interested in modifications that ex- 
plicitly break Lorentz invariance. From extrapolation of 
current experiments, we know that there exist stringent 
bounds on the most commonly expected types of Lorentz 
violations at the Planck scale (see, e.g., 0| and references 
therein). Nevertheless, even violations at much higher 
energy scales might still significantly affect black hole 
physics. 

In this paper we will investigate the effects that a su- 
perluminal modification of the dispersion relation would 
have on the Hawking radiation produced by collapsing 
configurations. Hawking's original prediction that black 
holes radiate thermally (T^, |Tl|, rested on the implicit 
assumption that the low-energy laws of physics, and in 
particular Lorentz invariance, are preserved up to arbi- 
trarily large scales, much higher than the Planck energy. 
The question of the robustness of Hawking's prediction 
to modifications of the transplanckian physics has been 
tackled principally by analyzing effective field theories 
such that at high energies a modification of the dispersion 
relation is incorporated (see however for a different 
take on the problem). Historically, attention has mainly 
been given to subluminal dis per sion relations. Impor- 
tant contributions such as [il,[ii,[il,[ii,[i3,[il[iltl3l, 
and more recently [2l|, [g^l seem to have settled the ro- 
bustness of Hawking radiation with respect to subluminal 
modifications and truncation of the frequency spectrum, 
although this conclusion still rests on certain assump- 
tions, usually related to the behaviour of the fields near 
the horizon. In any case, it is important to remember 
that this only solves part of the so-called 'transplanckian 
problem' with regard to Hawking radiation. Indeed, sub- 
luminal modifications gradually dampen the influence of 
ultra- high frequencies, and so they do not explore arbi- 
trarily large frequencies. So even assuming that it has 
been demonstrated that Hawking's result can be recov- 
ered without appeal to transplanckian frequencies, and 
that it is robust to subluminal modifications, the question 
remains whether it is also robust with respect to a (non- 
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dampened) modification of the transplanckian physics. 
Superluminal modifications differ conceptually from sub- 
luminal ones, in that they gradually magnify the influ- 
ence of ultra- high energies, and thereby offer an interest- 
ing test-case for the transplanckian robustness of Hawk- 
ing radiation. 

An additional motivation for our study comes from a 
possible connection with experiment. It is well-known 
(see, e.g., [l^) that Hawking radiation is not a purely 
gravitational effect, but a characteristic of quantum 
field theory in curved spacetimes with a horizon. Con- 
densed matter systems such as Bose- Einstein condensates 
(BECs) can, under certain approximations, be described 
by a relativistic quantum field theory. It has also been 
shown, again under certain idealizing approximations, 
that (acoustic) black hole horizons in BECs are dynam- 
ically stable [2J|^. Therefore, BECs are expected to be 
good candidates for a possible experimental verification 
of Hawking radiation ^6, 27, 28, 29,, 30, 31, M.M,M- 
But since the physics of BECs automatically leads to 
a superluminal dispersion relation at high energies, the 
question is again which kind of modifications are to be 
expected in the laboratory realization of a BEC black 
hole with respect to the standard Hawking picture. 

Of the above-mentioned works on the robustness of 
Hawking radiation, a few have also tried to address su- 
perluminal modifications [l^, [2l|, . Various problems 
make this case quite different from the subluminal one. 
These problems can be related to the fact that the hori- 
zon becomes frequency-dependent when modifying the 
dispersion relations. This is also true in the sublumi- 
nal case, though qualitatively in a very different way. 
With superluminal modifications, the horizon lies ever 
closer to the singularity for increasing frequencies, and 
asymptotically coincides with it. This causes the 'appar- 
ent' interior of the black hole (the interior of the zero- 
frequency horizon) to be exposed to the outside world. 
Since it seems unreasonable to impose a condition arbi- 
trarily close to the singularity as long as we do not have a 
solidly confirmed quantum theory of gravity, most of the 
approaches used for subluminal dispersion are invalid, or 
at least questionable, for superluminal dispersion. More- 
over, since it seems reasonable to expect that quantum 
effects will remove the general relativistic singularity, a 
critical frequency might appear above which no horizon 
would be experienced at all. 

In the analysis that follows we will try to avoid making 



Note, however, that this dynamical stability is quite fragile. For 
example, a minimal amount of backscattering in the interior of 
the black hole could be sufficient to create dynamical instabili- 
ties [24j . which in a realistic laboratory experiment might com- 
pletely obfuscate any quantum radiation related to an analogue 
Hawking process. This instability was first revealed in a mode 
analysis of a black hole-white hole configuration under a super- 
luminal dispersion relation |25|| . In this paper we will not worry 
about this issue and simply assume that the dynamics of the 
background is completely fixed. 



any further assumptions about the physics near the sin- 
gularity. For instance we will analyze the characteristics 
of the radiation at retarded times at which the singular- 
ity has not yet formed. Our approach will be based on 
a derivation of Hawking radiation through the relation 
between the asymptotically past ray trajectories and the 
asymptotically future ones in the case of a collapsing con- 
figuration. The only assumption about these asymptotic 
extremities are the standard ones, namely a Minkowski 
geometry in the asymptotic past, and flatness at spa- 
tial infinity also in the asymptotic future. The language 
we will use is related to the fluid analogy for black holes, 
which provides the intuitive picture of the spacetime vac- 
uum flowing into the black hole and getting shredded in 
its center. 



Our main results can be summarized as follows. Three 
crucial elements distinguish the late-time radiation with 
superluminal dispersion relations from standard relativis- 
tic ones. First, at any instant there will be a critical fre- 
quency above which no horizon has yet been experienced. 
This critical frequency will induce a finite limit in the 
modes contributing to the radiation, which will therefore 
have a lower intensity than in the standard case, even 
if the critical frequency is well above the Planck scale. 
Second, due to the effective frequency-dependence of the 
horizon, the surface gravity will also become frequency- 
dependent and the radiation will depend on the physics 
inside the black hole. Unless special conditions are im- 
posed on the profile to ensure that the surface gravity is 
nonetheless the same for all frequencies below the critical 
one, the radiation spectrum also undergoes a strong qual- 
itative modification. Depending on the relation between 
the critical frequency and the Planck scale, the radiation 
from high frequencies is no longer negligible compared 
to the low-frequency thermal part, but can even become 
dominant. This effect becomes more important with in- 
creasing critical frequency. Finally, a third effect is that 
the radiation will extinguish as time advances. 



The remainder of this article has the following struc- 
ture. In section |TT1 we will describe and motivate the 
classical geometry of our model and the concrete kinds 
of profiles that we are interested in. In section IIIIl we 
will briefly review how Hawking's standard result can be 
obtained in the formalism of our choice for the case of 
standard relativistic dispersion relations. This calcula- 
tion will be adapted in section lTVl to superluminal disper- 
sion relations, and we will analytically obtain a formula 
for the late-time radiation for this case. Then, in section 
IVl we will present graphics obtained from this analytic 
formula by numerical integration. These graphics will il- 
lustrate the results mentioned in the previous paragraph, 
which we will discuss in more depth and compare with 
other results in the recent literature in section IVll 
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II. THE CLASSICAL GEOMETRY 

We will study radiation effects in simple 1+1 dimen- 
sional collapsing configurations. It is well known that in 
the Hawking process, most of the radiation comes out in 
the s-wave sector. Hence, a spherically symmetric treat- 
ment, effectively 1-1-1 dimensional, suffices to capture the 
most relevant aspects of the process. We will work in 
Painleve-GuUstrand coordinates, where 

ds^ ^ ~[c^ ~v^{t,x)]dt'^ -2v{t,x)dtdx + dx^ , (2.1) 

and make the further simplification of taking a constant 
c (we will always write c explicitly to differentiate it 
from the frequency/ wave- number dependent Ck that will 
show up in the presence of a superluminal dispersion re- 
lation). In this manner all the information about the 
configuration is encoded in the single function v(t,x). 
The Painleve-GuUstrand coordinates have the advantage, 
compared to the Schwarzschild form, of being regular at 
the horizon. Moreover, they suggest a natural interpreta- 
tion in terms of the language of acoustic models. In such 
a model, c represents the speed of sound and v is the 
velocity of the fiuid (which corresponds to the velocity of 
an observer free- falling into the black hole). Throughout 
the paper we will consider the fiuid to be left-moving, 
f < 0, so that the outgoing particles of light/sound move 
towards the right. 

The most relevant aspects of the analysis presented 
hereafter depend only on the qualitative features of the 
fluid profile. However, to justify some specific calcula- 
tions later on it is helpful to use the following concrete 
profiles, see figs.[T]and[21 Consider a velocity profile v{x) 
such that v{x — > -l-oo) = 0, v{x = 0) = — c and further 
decreasing monotonically as x ^ —oo. In the fiuid image, 
the fiuid nearly stands still at large distances and acceler- 
ates inwards, with a sonic point or horizon at a; = 0. The 
further decrease in v can either be linear until a constant 
limiting value is almost achieved, as in fig. [1] or of the 
form 



v{x) 



2Af/c2 
+ 2Af/c2 



(2.2) 



again up to a constant limiting value, as in fig. [21 Let us 
recall that this second velocity profile corresponds to the 
Schwarzschild line element with M the central mass, as 
can be seen [35| by reparametrizing the time coordinate 



and using 



2M/c as radial coordinate. 



Up to here, we have described a stationary profile. 
Now, to incorporate the dynamics of the collapse, let 
us introduce a monotonically decreasing function £,{t), 
^{t —oo) +00, and define v{t,x) as 



v(t x) ^ - 



(2.3) 
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v(t) 
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FIG. 1: Velocity profile of a black hole with a linear slope 
of the velocity v{x), and hence a constant surface gravity 
K^/, from the classical horizon at a:: = down to some pre- 
defined limiting value. The auxiliary function ^(t) separates 
each instantaneous velocity profile v{t,x) into a dynamical 
region I: v{t,x) = v{t) — v{(,{t)) and a stationary region II: 
v{t, x) = v{x). 



to its gravitational (Schwarzschild) radius. The dynam- 
ical configuration that we obtain consists of a series of 
snapshots. In each snapshot, |w(t,x)| increases (i.e., the 
fiuid accelerates) from w = for a; ^ -|-oo, up to a 
point xq = £,(t), and then remains constant as x fur- 
ther decreases towards — oo. For consecutive snapshots, 
the point Xq = ^{t) moves leftwards, so v{t, x) covers an 
ever larger part of v(x). 

From the point of view of an outgoing (right-moving) 
particle of light/sound, the configuration is nicely split up 
into, first, a dynamical or i-dependent region, and then 




Imagining the collapse of a homogeneous star, the func- 
tion ^(i) represents the distance from the star surface 



FIG. 2: Velocity profile characteristic of a Schwarzschild black 
hole. The slope of v{x) increases from the classical horizon at 
X = leftwards up to a constant limiting value, leading to a 
surface gravity k.^i which increases with the frequency. 
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a stationary a;-dependent region, as defined in (|2.3p .^ 

The difference between tlie two types of profiles, figs.[T] 
andlH can best be explained in terms of the surface grav- 
ity. When modifying the dispersion relation, the horizon 
becomes a frequency-dependent concept: each frequency 
experiences a different horizon, as we will see explicitly 
in section IIVI In particular, for superluminal modifica- 
tions, the horizon forms later (i.e., at higher values of 
\v\, or more negative values of x) for increasing frequen- 
cies. The surface gravity will then also become frequency- 
dependent: 



dv 
dx 



(2.4) 



where xh,ui' is the frequency-dependent location of the 
horizon. This surface gravity will be seen in section IIVI 
and |V] to play a crucial role. This explains our choice of 
the two types of profiles: In our first profile, fig. [1] we 
have considered a linear velocity profile |w(a;)| from the 
classical or zero- frequency horizon {c — \v\) down to some 
given limiting value, |u|max; from where it stays constant. 
This maximum velocity defines a critical frequency 
as we will discuss in more detail in section IIVI Then 
the surface gravity k^i will be frequency-independent up 
to this critical frequency lo'^ after which it will rapidly 
vanish. In the second profile, see fig. [H we have taken a 
velocity profile typical for a Schwarzschild black hole, and 
therefore k^^' increases with lu' , again up to a predefined 
limiting value corresponding to the horizon for a critical 
frequency to'^. In the context of the fluid analogy, it seems 
obvious that some mechanism will avoid the formation 
of a singularity. But in any case, we will also take into 
account the limit for an infinite critical frequency, which 
corresponds to a velocity profile with |w|max ^ co- 



lli. STANDARD DISPERSION RELATIONS 

In this section we will briefly review a way of deriving 
Hawking's formula for the radiation of a black hole with 
standard (relativistic) dispersion relations of the form 
oj^ = c^k^. For the sake of simplicity, we will only con- 
sider a massless scalar field. We will summarize the main 
steps: from the Klein-Gordon inner product and the def- 
inition of the Bogoliubov f3 coefficients, over the relation 
between past and future null coordinates, to the black- 
body radiation in the wave-packet formulation. Our aim 
is to present the key points of the procedure in such a way 
that they can easily be adapted to superluminally modi- 
fied dispersion relations — the subject of the next sections. 



In a physically realistic model, the apparent kink in the profile 
where the transition between both regions takes place will of 
course be smoothed out. 



A. Inner product 

The d'Alembertian or wave equation for a massless 
scalar field in 1-1-1 dimensions for the metric (|2.ip with 
constant c can be written as 



{dt + dxv){dt + vdx) 



c'd: 



(3.1) 



This conformally invariant theory is equivalent to the 
dimensionally- reduced 3-1-1 spherically-symmetric theory 
if one neglects the backscattering due to the angular- 
momentum potential barrier (responsible for the so- 
called grey-body factors). 

In the space of solutions of this wave equation, we can 
define the Klein-Gordon pseudo-scalar product 



d^^ ipi 



(3.2) 



which is independent of the choice of the spatial slice E. 
For a t =constant slice, and in particular for t -l-oo, 
this becomes 

(<y5i,<i52) = \ dx \px{dt + vd^)kp2- ^*2{dt^vdx)'f-\\ . 



(3.3) 

We can define future null coordinates u(i, x) and w(t^ x), 
such that, when t,x +oo, 

u{t, x) ^ t — x/c , w{t, x) ^ t + x/c . (3.4) 

Writing the inner product in terms of these null coor- 
dinates, gives, for the limit t — > -foo. 



+ 00 



du [ipidu^p*2 - ^Idu^pi 



\w—'\-oo 



+ 00 



(3.5) 



Note that, in the derivation of this formula, we have only 
made use of coordinate transformations, without any ap- 
peal to their null character or to geometrical tools such 
as conformal diagrams or the deformation of Cauchy sur- 
faces. This is important because in the case of modified 
dispersion relations, such geometrical concepts become 
problematic and actually can only be maintained in the 
context of a rainbow geometry — if at all. 

Similarly we could have used past null coordinates 
([/, W), which obey, when t — > —00, 

U{t,x)~^t~x/c , W{t,x) ^t + x/c , (3.6) 

to calculate the inner product in the asymptotic past. 

B. Bogoliubov /3 coefficients 

The right-moving positive-energy solutions associated 
with the asymptotic past and with the asymptotic future, 
normalized in the Dirac-delta sense, can be expressed as 



1 



1 



(3.7) 
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respectively, where we use primes to indicate asymptotic 
past values. 

The relevant Bogoliubov f3 coefficients encoding the 
production of radiation are defined as Pi^^j' = (V'L'' V'^)- 
The mode mixing relevant for the Hawking process occurs 
in the right-moving sector. Therefore, we only need to 
calculate the first term in the scalar product (13. Sp . Then, 
plugging p.7p into the definition of f3 and integrating by 
parts we obtain the simple expression 



B ,- — I- 



(3.8) 



so that all the information about the produced radiation 
is contained in the relation U = U (u) = U{u, w +00). 



Then, from the expression 

/3u,j,ur,i^' = J '^^ Pu,u,'Pi^,.ui{'^) , and assuming that 
the wave packets are sufficiently narrow (Aw ^ ujj), we 
obtain 



cAu! sin^(z — zi) 



1 



(3.13) 



2TrLu'K {z-ziY exp(2l£^)-l 

where we have defined 

cAuj , Auj 
z=— — inuj Ao , zi^—iui-ui). (3.14) 

Finally, integration in oj' gives the number of particles 
with frequency Wj detected at time ui by an asymptotic 
observer: 



C. Relation U{u) 

For a standard relativistic dispersion relation, it is well 
known that the relation between U and u for a configu- 
ration that forms a horizon (in our case at a; = 0) can 
be expressed at late times as f7 = Uh — ^le"""/^ (note 
that we use a subscript H for all quantities associated 
with the horizon), where Ur^ A and the surface gravity 



dx 



(3.9) 



x=0 



are constants. 

We can define a threshold time ui at which an asymp- 
totic observer will start to detect thermal radiation from 
the black hole. This retarded time corresponds to the 
moment at which the function U (u) enters the exponen- 
tial regime. We can then rewrite the previous expression, 
valid for w > w/, as 



U = Uh- Ane 



— k{u — Uj ) /c 



(3.10) 



Plugging this relation into (j3.8p and integrating in u gives 

1 C I CLjJ , 

Pujuj' = — W — - exp[-«a; [/ffj exp[-i — ln(w Aq)] 

Zt: \ LU K K 

X exp(— iwM/) exp( — - — )T{icuj/K) . (3.11) 



D. Wave packet formulation 

In order to obtain physically sensible results, it is a 
good precautionary measure to replace the monochro- 
matic rays described until now by wave packets (see 
e.g. [ll| or the discussion in ^36^). Positive-energy wave 
packets can be defined as 







otherwise. 



< (3.12) 



where ui = uo + 2ttI/Auj with uq an overall reference and 
I an integer phase parameter. The central frequencies of 
the wave packets are LUj = jAuj, with Auj their width. 



N,, 



UJj,Ui 



exp {2Tr cuj j / k) — 1 



(3.15) 
(3.16) 



which reproduces Hawking's formula (in the absence of 
backscattering) and corresponds to a Planckian spectrum 
with temperature Th = k/ {2-kc). 



IV. SUPERLUMINALLY MODIFIED 
DISPERSION RELATIONS 

In this section we will indicate how the late-time radi- 
ation originating from the formation of a black hole from 
collapse can be calculated in the case of superluminal 
dispersion relations. 

We introduce superluminally modified dispersion rela- 
tions by adding a quartic term to the wave equation: 



{dt + d^v)(dt + vd^)^^^ 



(4.1) 



where kp (the 'Planck scale') is the scale at which non- 
relativistic deviations in the associated dispersion rela- 
tion 



(4.2) 



become significant. We use this relation for concreteness 
and because it is the one that shows up in Bose-Einstein 
condensates, where kp = 2/^, with ^ the healing length 
of the condensate. However, qualitatively, our results will 
not depend on the specific form of the deviations from the 
relativistic dispersion relation but on their superluminal 
character. 

This dispersion relation leads to a modification in both 
the phase velocity Vph and the group velocity Vg. For a 
right-moving wave, we have 



Vph = T 



Ck,ph 



= dfc " 



(4.3) 
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FIG. 3: Behaviour of the effective phase Ck,ph and group Ck,g 
speeds of light/sound with respect to the wave number k. 
Due to the non-relativistic dispersion relation, the effective 
velocities become 'superluminal' for k > kp (where kp is the 
Planck scale). 



where we have introduced the effective fc-dependent 
phase and group speeds of Ught/sound 



/ fc2 
Ck,ph = cW 1 + ^ 



1 



Ck,. 



(4.4) 



respectively (see fig. [3]). Both Ck,ph and Ck,g become 
larger than c (i.e., become 'superluminal') as k increases 
above kp. 

At first sight, it seems obvious that the ray equation 
should be defined in terms of the group velocity Vg . Nev- 
ertheless, the question of whether the velocity relevant 
for Hawking radiation is the _group or the phase velocity 
seems to be tricky [g^, [s^, Isl]. For example, (|4.9p be- 
low suggests that the phase velocity might be relevant. 
We limit ourselves to remark that Ck,g and Ck,ph show 
the same qualitative behaviour. Hence our results are 
independent of this question and we will simply write 
Cfc (or Ck{uj'} when wishing to emphasize the frequency- 
dependence). Then, there will be a horizon, which be- 
comes frequency-dependent, when Ck + v = 0, irrespec- 
tively of whether Ck,g or Ck^ph is used for Ck- More- 
over, since becomes arbitrarily high for increasing wave 
number A:, given a certain |w|rnax at a particular instant 
of time, there will be a critical such that waves with 
an initial frequency oj' > uj'^ do not experience a hori- 
zon at all. The only exception to this rule occurs when 
the velocity profile ends in a singularity v —oo, which 
implies — > oo. 

Our aim is to calculate the black hole radiation with su- 
per luminally modified dispersion relations. We will now 
repeat the main steps of section Hill and point out where 
and how these modifications must be taken into account. 



A. Generalization of inner product 

The essential point with regard to the pseudo-scalar 
product l|3.3p is the following. Its explicit form for 
t =constant is not changed by the presence of the 9^ 
term in the wave equation (see also the discussion in [25j). 



Indeed, it is still a well-defined inner product, and in par- 
ticular it is still a conserved quantity, since 



9t(<^l,</52) = dx ipi dl , 



(4.5) 



which can be seen, by repeated integration by parts, to 
vanish under the usual assumption that the fields die off 
asymptotically. Note that the modification of the dis- 
persion relation singles out a preferred time frame: the 
'laboratory' time t. Changing to another time t will in 
general lead to a mixing between t and x, and hence the 
simple relations given here would no longer be valid. 

Using the preferred time i, and making exactly the 
same change of coordinates as in the case of standard 
dispersion relations, we can again transform the inner 
product, evaluated at t -f cx), into the expression (j3.5p 
in terms of u and w. As in the standard case, only the 
first term is relevant for the Hawking process: 



"2 



+ 00 



\'W — ArOO 



(4.6) 



Note that we are now using u and w merely as a per- 
fectly good set of (auxiliary) coordinates, in order to cast 
the inner product into a useful form. However, we cannot 
associate them any null character such as they had in the 
case of standard dispersion relations. 



B. Rainbow null coordinates 

Let us define some sets of spacetime functions that will 
prove to be useful in what follows. Given a fluid profile 
of the type described in section |TT1 we can integrate the 
ray equation 



da; 



Ck(uj'){i,x) -f v(t,x) , 



(4.7) 



starting from the past left infinity towards the right. 
The ray initially has a frequency w' and an associated 
initial wave number k' , from which we can deduce the 
value of Ck' ■ In the left region, where the velocity profile 
is dynamic but position-independent (v{x,t) = v{^(t)), 
with £^(t) the auxiliary function introduced in section ITT)) . 
k = k' can be considered as fixed while the frequency 
changes (this is what happens in a mode solution of equa- 
tion (|4.ip in this region). Then, we can define the func- 
tion 



If k' 

U^'{t,x) = — / oj{t)dt-—a 

UJ' J UJ' 



(4.8) 



where ijj{t) is the instantaneous frequency of the particle 
at each time t, defined through the dispersion relation 
and such that id{t — oo) = w' . 

When the ray reaches the kink it passes into a station- 
ary region in which the velocity profile only depends on 
the position {v{t,x) = w(a;)). At the kink, the ray still 
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has the initial wave number k' and a frequency lj. In 
its propagation towards the right this frequency now re- 
mains fixed while its wave number becomes a decreasing 
function k{x) of the position, such that fc('^(io)) — k', 
with to the moment at which the kink is crossed. The 
final frequency of the ray will then simply be lu and its 
final wave number k = limx^+ook{x). In this region, 
then, the function U^i {t, x) can be expressed as 



k{x) 



dx . 



(4.9) 



Note here that uj and k{x) both depend on the initial lo' . 

The same can be done by integrating the ray equa- 
tion starting from the future. In this way we can define 
functions. (The same procedure can be used to de- 
fine Wuj' and Wui). It is worth noting that lAuj' and Uuj 
are not null coordinates in the usual geometric sense, 
since they are frequency-dependent, but they neverthe- 
less share many properties with null coordinates. 



C. Bogoliubov 13 coefficient 



The calculation of the inner product (|4.6p involves the 
limit w +00 (equivalently, Wuj — > -foo). So, the gen- 
eral form 

= u^{u,w) (4-10) 

can be simplified to u^^ = u^{u) = Uuj{u,w +oo). 
We can therefore change variables in the inner product 
from u to Ui^. The combination of the derivative and 
the integral means that the form of (|4.6p is preserved 
in the new integration variable u^j- Then we can write 
the Bogoliubov (3 coefficients relevant for the Hawking 
process as 



(4.11) 

Now, assuming that the profiles vary slowly (in scales 
much larger than the Planck distance) , the right-moving 
positive-energy modes associated with past and future 
infinity can be approximated by the following simple ex- 
pressions: 



1 



^ — iuj' lA, , 



1 



so that the Bogoliubov coefficients read 



1 ho 



o'U^f iu^) —iujn 



(4.12) 



(4.13) 



(4.14) 



In analogy with the standard case, all information 
about the radiation is seen to be encoded in the relation 



l^uj'{ucj). In this expression uj' is the initial frequency of 
a ray at the past left infinity and uj = w(w') its final 
frequency when reaching the future right infinity. 

The previous approximation for the Bogoliubov coeffi- 
cients amounts to considering profiles that vary slowly 
both with X and t. This is equivalent to considering 
large black holes. In general the quartic term in the wave 
equation, or equivalently the quartic modification of the 
dispersion relation, introduce a new source of backscat- 
tering, on top of the usual angular-momentum potential 
barrier which we have already neglected. In our approxi- 
mation, this additional backscattering (beyond the stan- 
dard grey-body factors) has been neglected. For large 
black holes this contribution will in any case be very small 
as has been observed in numerical simulations f25| . In ad- 
dition we are also neglecting any reflection caused by the 
kink. However, let us remark that given an approxima- 
tive scheme for calculating Ul^'{u^) for general profiles 
v{t,x), the same lA^i{u^) obtained from a profile with a 
kink would be obtainable from one (or several) specific 
v{t,x), this time perfectly smooth and thus causing no 
further backscattering. Our results, which rely only on 
the specific form of the relation lA^i(u^), are therefore 
valid beyond the specific configurations with a kink pre- 
sented in this paper. 



D. Relation U^i{uuj) 

Our next task is to calculate the uniparametric family 
of functions 14^^' and the relation between U^^i and for 
different configurations. As explained in section IIVBI 
the relation lAuj'{ui^) is obtained by integrating the ray 
equation dx/dt = Ck + v using the profiles discussed in 
section [TTl These can be described by means of station- 
ary profiles v(x) and an auxiliary function ^(t), see the 
definition of v{t^x) in (|2.3p . We will use a straightfor- 
ward extension of the procedure established in [39| (see 
also [40j for a summary) for a relativistic dispersion re- 
lation. Care must be taken, however, with the quantities 
that depend on the frequency lo' . In particular: 

• We will denote by xh,uj' and tn^^j' the position and 
the time at which the horizon associated with a 
particular initial frequency uj' is formed. 

• The surface gravity k^/, defined in (|2.4p . allows to 
write, for all uj' < uj'^ and for x close to xh : 



1 



V{x) W -Cfe + -K^'{X - XH,uj') 

c 

up to higher-order terms in x — xh,u 



(4.15) 



• The linearization of ^(i) near tu^u)' also requires the 
introduction of an cj'-dependent parameter X^^r. 



S_{t) « XH,uj' - \oj'{t - tH,uj') 

again up to higher-order terms. 



(4.16) 
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Let us consider laboratory times t > tH.o = tH,u'=o 
such that ^ {t) has already crossed the classical horizon at 
X — Q: ^{t > tn.o) < 0- As we explained earlier we can 
define a critical frequency uj'^ as the minimum initial fre- 
quency such that Ck + v{^{t)) > for all uj' > cu'^, or, 
in other words, the minimum frequency which at that 
particular time has not experienced any horizon yet. 

In the dynamical part of the profile, where v = v{^{t)), 
for rays crossing the kink just before the formation of the 
horizon, integration of the ray equation leads to 



Ck' 



(to, 



tH,L 



(4.17) 



for small values of |io,w' — i_f/,w'| (see [39'] for details of 
the calculation in the case of a relativistic dispersion re- 
lation), where fo.tj' is the (frequency-dependent) time at 
which the kink separating the dynamical and the station- 
ary regions is crossed, Ck' = lim(^_oo Ck{uii) is the speed 
of light in the asymptotic past, and 



U 



^ tH.u.' - + — f"'^' mmt (4.18) 

Cfc' Ck' J_oo 



is the ray constituting the horizon associated with the 
frequency lu' . 

In the stationary part of the profile, where v = v{x), 
we obtain, for io,w' — > tn^i^', again using a method and 
notation based on (39i] . 



e(to,^') 



■ lnK(io, 



XH.i 



CkJ Ku 

(4.19) 

where Ckj = limt^+oo Ck(uj') is the speed of light in the 
asymptotic future, and Ci a bulk integration constant. 
So 



(4.20) 



Both regions are connected by combining (|4.17p and 
(13201) • Making use of (gUl), this gives 



(4.21) 



This relation is valid for frequencies for which a hori- 
zon is experienced, i.e. for ld' < oj'^ and times > 
■"/.oi'i where uj ^^i is the threshold time defined in sec- 
tion nil CI (which, unsurprisingly, has become frequency- 
dependent). Again, as in the case of standard dispersion 
relations — section IIII CI — we can write 



U,., = U 



-Aoe 



,)/c 



(4.22) 



where ui^^' is essentially uj^^^i (with possible higher-order 
corrections). Assuming that the collapse takes place 
rapidly, it is a good approximation to replace ui^uj' by 
^i,uj'^- Actually, as we will see shortly, this is a conserva- 
tive estimate, in the sense that it slightly underestimates 
the superluminal correction to the radiation spectrum. 



E. Wave packet formulation 



The relation Uuj'(u^) we are considering interpolates 
between a linear behaviour at early times and an expo- 
nential behaviour at late times. It has the same form as 
the relation U{u) for standard dispersion relations, and 
so we can continue following the steps of the standard 
case. In particular, the equivalent of p. lip is obtained 
by integrating out u^- Note that we use the subscript uj 
to emphasize that the u^^ are not the null coordinates of 
the standard case, but this should not be interpreted as 
an explicit function of to and so does not complicate the 
integration steps in u^^ and w. However, when integrat- 
ing |/3tjj,u!;w'P in uj' to obtain N^-^m, see section ITlI Dl 
we must carefully consider the frequency-dependence of 
the relevant terms, i.e., ofUH,uj', and k^/. The term 
carrying L(h,uj' is moduloed away in p.l3p . and we have 
replaced ui.uj' by ui.u'^, so the only relevant frequency- 
dependent factor that we are left with is the surface grav- 
ity . 

Moreover, because of the critical frequency to'^ in the 
horizon formation process, a finite upper boundary will 
also be induced in the integral. Indeed, frequencies 
Lu' > u'^ do not contribute to the radiation at all, since 
they do not experience a horizon. This is a delicate but 
crucial point. It was already observed long ago by Jacob- 
son [14] that trying to solve the transplanckian problem 
naively by imposing a cut-off frequency would seemingly 
extinguish Hawking radiation on a relatively short time 
scale. In our case, however, this cut-off is not imposed ad 
hoc, but appears explicitly because of the superluminal 
character of the system at high frequencies. Moreover, 
the critical frequency, and hence the upper boundary 
induced in the integral, depend directly on the physics 
inside the horizon. Indeed, given a certain velocity pro- 
file, and in particular its behaviour near the center of the 
black hole, the critical frequency can be calculated by set- 
ting Ck = \v\ in ea. (|4.4p and extracting the corresponding 
critical frequency from the dispersion relation (|4.2p . We 
will see this effect graphically in section IVl 



In analogy with eqs. p.l3p and (|3.15p . we now obtain 
the number of particles detected for each frequency uij as 



■ dc^' 1/3, 



cAcu 
"2^" 



duj' 1 



UJj ,Ui 

sm 



(4.23) 



W 



exp(-;-T^) - 1 



where now z = Incj'An, with kq = k^'=o (which 

ZAvQ ^^^^ 

corresponds to the standard k, of p.9p ). and 
zi,uj' = '^^{ui — uj^^i^). Changing the integration 
variable from uj' to z, we finally obtain the central ex- 
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pression in our analysis: 



N, 



UJj ,ui 



/ ZTTCOJj \ 1 



(4.24) 

Note that, as we indicated earlier, the use of u/.c^' in the 
definition of zi^^j' is a conservative stance. Indeed, strictly 
speaking, we should write zi^^^i (x {ui — ui.uj>). For a fixed 
Zc, a smaller value of z;,cj', and hence a larger value of 
ui^uj'i rneans that a larger part of the central peak of the 
integrand will be integrated over. Since we are replacing 
ui,ui> by the upper bound u/^^^, we are overestimating 
the resulting radiation (i.e., underestimating the modifi- 
cation with respect to the standard Hawking radiation). 

The expression (|4.24p brings out the two crucial factors 
mentioned earlier, and a third, corollary one. 

• First, it shows the dependence of the total radiation 
on the critical frequency uj'^ (through the integra- 
tion boundary Zc induced by it), as discussed just 
before obtaining formula (|4.23[) . 

• Second, it shows the importance of the frequency- 
dependent Kuj' (to be compared with the fixed k of 
the standard case). Given a concrete profile v{t, x), 
the frequency-dependence of Kai' can be derived ex- 
plicitly, as we will discuss next. 

• Finally, as a corollary of the first point, it shows 
that, as ui (and hence zi^^i) increases, a smaller 
part of the central peak of the integrand will be 
integrated over, so the radiation will die off as it; 
advances. 
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FIG. 4: Surface gravity k^^/ with respect to the (logarithm 
of the) frequency for a Schwarzschild-type black hole as in 
fig. O with the Planck scale kp — 10^^ and the critical 
frequency lo'^ — 13x 10'^^ t~^ (where t represents an arbitrary 
time unit.) 



Given a concrete profile, Xh.u' can then be obtained and 
Ki^' calculated explicitly. 

For the profiles of the first type, see fig. [1] the result 
is trivial. Indeed, |w(a;)| increases linearly between the 
horizon corresponding to w' = and the one for w^, so 
we obtain a constant k^j' = for all lu' < u>'^. 

For a Schwarzschild profile as in fig. ^ on the other 
hand, 



vlx) 



I cV2^o 
X + c^/2ko ' 



(4.27) 



so we obtain 



F. Surface gravity 

A careful analysis has shown that, so far, most formu- 
las for the standard case could be adapted to superlumi- 
nal dispersion relations by replacing the relevant magni- 
tudes with their frequency-dependent counterparts. For 
example, P^ju' was obtained by replacing U and u by L(i^> 
and Uuj, and in particular Uh, ui and k by Uh,uj', uj^^r 
and , respectively. Given a concrete profile, we can ex- 
plicitly deduce the relation between Ht^i and to' , with the 
Planck scale fcp as a parameter, as follows. The horizon 
for a particular initial frequency uj' is formed when 



1 + % 
k p 



(4.25) 



i.e., when Vph = Ck + v = 0, where we have used the 
phase velocity for concreteness. But again, qualitatively, 
our results would be similar if taking the group velocity 



Vg instead of Vph- 



Taking into account that kn = fc', the dispersion rela- 
tion in the asymptotic past can be written as 



- 1 



(4.26) 




(4.28) 



Note that Ki^i remains nearly constant until frequencies 
of the order of magnitude of the Planck scale are reached, 
and then starts to increase rapidly, see fig. [H As we will 
see graphically in the next section, this can have an im- 
portant qualitative influence on the radiation spectrum. 

We now have all the tools necessary to compute and 
plot (|4.24p for different parameters of the profiles de- 
scribed in section HIl 



V. GRAPHICAL RESULTS AND DISCUSSION 



We have calculated eq. (|4.24p numerically using the 
Gauss- Chebyshev quadrature method. The results are 
plotted in figs. [SHU and perfectly illustrate the three 
important factors that we deduced theoretically in sec- 
tion IIVEI Note that in all the figures we have plotted 
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FIG. 5: Influence of the critical frequency on the radiation 
spectrum for a black hole with velocity profile such that 
the surface gravity is constant, as in fig. [T] and different 
values of the critical frequency. We have chosen c = 1, 
ui = ui = 0, Uh — a = It, k.0 = 12 and the 
Planck scale kp — 10^^ t~^ (where t represents an ar- 
bitrary time unit), which amounts to considering a solar- 
mass black hole. From top to bottom we have plotted 
Lj', = 103°°°,10l°^103^1,10-='^10-™,10-3°°. Note that 
the standard Hawking spectrum coincides perfectly with the 
upper curve, which effectively corresponds to the absence of 
a critical frequency. 



E — uj^ X N against uj to make visual comparison with 
the usual thermal spectra in 3-t-l dimensions easier. 

Fig. [5] shows the influence of the critical frequency 
for a profile with constant k^i , as in fig.[Tl and ui — ui^^j' 
(i.e., immediately after the horizons have formed for all 
Lo' < uj'^. The constant k^i' guarantees that the form 
of the thermal spectrum is preserved. However, the in- 
tensity of the radiation decreases with decreasing crit- 
ical frequency. Actually, only for extremely high criti- 
cal frequencies is the original Hawking spectrum recov- 
ered. For a critical frequency still well above the Planck 
scale, the decrease can be significant. For example, for 
Lo'^ — lO^^c/cp, the peak intensity decreases by nearly 
30%, while for lo'^ = ckp, the decrease is approximately 
40%. At the other end, note that extremely low critical 
frequencies still leave a significant amount of radiation. 
For example, for u)'^ = 10~^^^c/cp, stiU 20% of the orig- 
inal peak radiation is obtained. So Hawking radiation 
receives significant contributions from an extremely wide 
range of frequencies. 

Fig. [S] shows the radiation spectrum for a 
Schwarzschild-like profile and hence increasing k^i, 
as in fig. [21 in particular for critical frequencies lj'^ close 



FIG. 6: Influence of a frequency-dependent surface grav- 
ity Hi^i on the radiation spectrum for a black hole 
with Schwarzschild-type velocity profile (surface gravity in- 
creases with frequency), as in fig. [51 and different val- 
ues of the critical frequency around the Planck scale kp: 
uj'c/ ckp — 13,10,7,4,0.1 (from top to bottom). The standard 
Hawking spectrum is depicted in dashed-dotted line for com- 
parison. Numerical values of kq, kp etc as in fig. [5l 

to the Planck scale. On top of the general decrease of the 
standard thermal part of the spectrum (approximately 
40%, as in the previous case) due to the finite integration 
boundary induced by the critical frequency, the fact that 
the surface gravity is now frequency-dependent leads 
to an important qualitative change of the spectrum. 
The high-frequency tail of the spectrum is totally trans- 
formed. Actually, if the critical frequency is sufficiently 
higher than the Planck scale, the dominant source of 
radiation lies in the high-frequency region. Note that 
this effect is truly a consequence of the modification of 
the physics for frequencies above the Planck scale. This 
can be appreciated by noticing that, at u'^ = O.lcfcp, 
the whole tail-modifying effect has disappeared and 
the usual thermal form of the radiation spectrum is 
recovered (although still with the quantitative decrease 
described above). 

Finally, fig.[71shows the influence of the measuring time 
ui (measured with respect to uj^i^i ) for a profile of the 
second type (increasing k^'). It is clearly seen that the 
radiation dies off with time, and actually dies off rather 
fast. For an actual solar-mass black hole and an oj'^ of 
the order of the Planck scale, the radiation would last 
only a few tens of milliseconds. Note that this effect is a 
corollary of the existence of a critical frequency w^, since 
in its absence, the integration boundaries in (|4.24p would 
be infinite, and so the integral would be insensitive to a 
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FIG. 7: Influence of the measuring time Ui on the radiation 
spectrum, for a Schwarzschild-type velocity profile, see fig. [5] 
and uj'c/ckp = 13, compare with fig. ID Different values of ui 
(from top to bottom): it; — 0, 22, 35, 50, 500 t with t denoting 
an arbitrary time unit. The standard Hawking spectrum is 
depicted in dashed-dotted line for comparison. Numerical 
values of no, kp etc as in fig. [S] 

change u u + Au. 

The combined effect of fig. IH] and fig. [7] leads to the 
following qualitative picture for the further collapse to- 
wards a singularity once the initial or classical horizon 
has formed. As the interior gradually uncovers a larger 
portion of the Schwarzschild geometry, two competing 
processes will take place. On the one hand, the spectrum 
acquires ever larger contributions associated with higher 
and higher temperatures. On the other hand, the overall 
magnitude of the spectrum is damped with time. The 
question of which process dominates would depend on 
the fine details of the dynamics of the collapse, and might 
be further complicated by backreaction effects, which we 
have not considered in our analysis. 

VI. SUMMARY AND CONCLUSIONS 

We have discussed the Hawking radiation for a collaps- 
ing configuration with superluminal dispersion relations. 
Modifications of the dispersion relation cause the hori- 
zon, and various associated quantities such as the surface 
gravity k, to become frequency-dependent. In particular, 
a critical frequency uj'^ naturally appears such that fre- 
quencies higher than uj'^ do not experience a horizon at 
all. More generally, it also means that the standard ge- 
ometric concepts traditionally used to study black holes 
must be handled with care. Nevertheless, through a de- 



tailed analysis, we have seen that the equations related to 
the late-time radiation can be adapted quite straightfor- 
wardly from standard (relativistic) to superluminal dis- 
persion relations. 

We analytically derived an approximate equation for 
the particle production at late times (|4.24p with super- 
luminal modifications of the dispersion relations. This 
equation clearly showed that important modifications in 
the late-time radiation should be expected, first, due to 
the existence of the critical frequency tu'^ and the finite 
upper boundary it induces in the integration, and second, 
due to the frequency-dependence of k^'. We integrated 
(|4.24p numerically, and plotted the resulting spectrum, 
thereby confirming these expectations. 

We have seen that the standard Hawking spectrum 
is recovered only in a very particular case: in the limit 
when the critical frequency goes to infinity (i.e., when 
the profile for the velocity |w| goes to a singularity) and 
moreover the surface gravity Ki^> is constant (linear ve- 
locity profiles of the type of fig. [J). For lower critical 
frequencies, as long as k remains constant, the thermal 
form of the Hawking spectrum is maintained, but the 
intensity decreases rapidly with decreasing lj^. A non- 
negligible radiation persists, however, even for extremely 
low critical frequencies. This clearly establishes the first 
point, namely the importance of the critical frequency w^, 
which illustrates the statement made in the introduction 
that with superluminal dispersion relations, the interior 
of the black hole is probed, which can significantly affect 
the Hawking process. 

For a Schwarzschild-like profile, the velocity profile 
generally leads to a frequency-dependent surface grav- 
ity Koj' , increasing with the frequency lu' . This means 
that the standard Hawking result cannot be recovered for 
a Schwarzschild black hole with superluminal dispersion 
relations, even when the velocity profile has a singularity. 
Actually, when the velocity profile has a limiting value, 
the same quantitative decrease of the Hawking part of 
the spectrum as before shows up. Moreover, if the criti- 
cal frequency oj'^ is above the Planck scale, a drastic qual- 
itative change of the radiation spectrum takes place, and 
for sufficiently high uj'^ (roughly a few times the Planck 
scale) the high-frequency part of the spectrum even be- 
comes dominant. This shows the importance of the sur- 
face gravity, which again illustrates the role played by 
the interior of the black hole. 

Finally, we have also seen that, as a corollary of the 
existence of a critical frequency and of the finite upper 
boundary induced by it in the thermal response function, 
the radiation spectrum dies off as time advances. 

A few observations might be useful to connect our work 
with existing results on Hawking radiation and its sensi- 
tivity to modified dispersion relations. The most general 
observation is that the 'robustness' of Hawking radia- 
tion, which is often considered to be a well-established 
result, actually depends crucially on a series of assump- 
tions. These assumptions might be reasonable in the case 
of subluminal dispersion relations. But for the case of su- 
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perluminal dispersion relations, as we have shown explic- 
itly, the assumptions needed to reproduce the standard 
Hawking result depend on the physics inside the (zero- 
frequency) horizon, and moreover in a way which is not 
compatible with the Schwarzschild geometry. In partic- 
ular, in [20l |. it was shown for a stationary scenario that 
Hawking radiation is robust with respect to superluminal 
modifications of the dispersion relations, provided that 
positive free fall frequency modes were in their ground 
state just before crossing the horizon. However, it was 
also admitted that it is not clear whether this is the phys- 
ically correct quantum state condition. In [2l| . three ex- 
plicit assumptions were given for the previous condition 
to hold: freely falling frame, ground state and adiabatic 
evolution. In the same article, it is similarly pointed 
out that these assumptions might fail for a superluminal 
dispersion relation, since (as we also mentioned in the 
introduction) for high-frequency modes this implies that 
one makes an assumption about the physics at the singu- 
larity. Rather than imposing any conditions on the state 
near the horizon, and hence ultimately at the singular- 
ity, we have fixed the initial geometry to be Minkowskian 
and the quantum field to be in the natural associated 
vacuum state, and evolved this into a black hole config- 
uration. Since our results are seemingly in contradiction 
with the ones obtained in po| and [2l|, it is worth ex- 
amining in which sense the conditions stated by those 
authors are violated in our approach. If one considers a 
collapse scenario, for example of a BEC, in a laboratory 
setting, then it is quite natural to assume that the "freely 
falling frame" condition will be violated. Indeed, the su- 
perluminal modification is in this case associated with 
the existence of a privileged external reference frame: 
the black hole rest frame or lab frame (as we noted in 
sec . Iiy A|) . and not the freely falling frame, as assumed 
in [2l|. Note that this violation of the free-fall frame 
condition automatically implies a violation of the ground 
state condition in the sense in which this condition is 
formulated in [21], namely as the ground state "with re- 
spect to the freely falling frame" . Indeed, our results sug- 
gest the following picture. The low-energy modes experi- 
ence the classical geometry, are therefore dragged along 
in the free-falling frame and so at the horizon they occupy 
the vacuum state associated with this free-falling frame, 
namely the Unruh vacuum. Hence they contribute to 
the black hole radiation in the traditional Hawking way. 
However, the ultra-high-frequency modes (above the crit- 
ical frequency) do not see the horizon. Hence they do 
not couple to the classical geometry of the collapse, but 
rather remain connected to the external or laboratory 
frame, and therefore pass through the black hole (nearly) 
undisturbed. The ground state of these high-frequency 
modes is then not the vacuum associated with the freely 
falling frame, but the Boulware vacuum associated with 
the original Minkowski geometry, or in other words, with 
the stationary reference frame of the lab. So these fre- 
quencies above the critical frequency do not contribute 
to the thermal output spectrum. The overall radiation 



is a convolution of the contributions from all the differ- 
ent initial frequencies, where the surface gravity for each 
frequency can be interpreted as leading to an effective 
weighting factor. Depending on the internal physics of 
the black hole, this leads to either a spectrum of the 
traditional Hawking form but with a reduced intensity 
(in the case of a constant surface gravity) , or a modified 
spectrum where the high-frequency contributions dom- 
inate (as for a Schwarzschild profile, provided that the 
critical frequency lies sufficiently above the Planck scale) . 
In any case, what our analysis shows with respect to the 
situation treated in [l^ and [2l| is the following. The 
robustness of Hawking radiation discussed by those au- 
thors is the robustness with respect to a superluminal 
modification under the precise conditions occurring in 
the standard case, in which all modes occupy the ground 
state associated with the free-falling frame just before 
crossing the horizon. This assumption implies that all 
modes, regardless of their frequency, occupy the Unruh 
vacuum. Our analysis shows precisely that the super- 
luminality can lead to a spontaneous breaking of these 
conditions. The natural evolution of a collapsing con- 
figuration is such that the frequencies above the critical 
frequency are not in the Unruh vacuum associated 
with the freely falling frame, but remain in the Boulware 
vacuum associated with the initial asymptotic Minkowski 
condition. It should then come as no surprise that this 
can indeed have an important impact on the resulting 
radiation spectrum. 



We should mention that in a recent paper by 
Schiitzhold and Unruh [22] based on an alternative ap- 
proximation method, a new potential problem with su- 
perluminal dispersion relations is highlighted. This 
"ultra-violet catastrophe" is associated with a higher- 
than-linear growth of the group velocity at large wave 
numbers, and should therefore not occur for the disper- 
sion relation used in our work. However, we have found 
that a Schwarzschild interior geometry could give place 
to the same sort of phenomenology. 



A final observation concerns the recent article [33[, in 
which the authors numerically simulate the formation of 
an acoustic horizon in a BEC, and analyze the creation 
of correlated pairs of phonons through the so-called trun- 
cated Wigner method. Our findings for the case of a 
constant surface gravity seem to be in qualitative dis- 
agreement with the discussion presented in ^33j ■ since 
the authors of this paper assert having observed a sta- 
tionary Hawking flux, while our analysis concludes that 
the radiation should quickly fade off. The source of this 
apparent discrepancy might reside in the fact that the 
correlation function which they study is normalized, and 
therefore probes the form of the spectrum (which we also 
find stationary with time), but does not provide details 
about the net amount of particle production. 
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